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SCALING LIMIT FOR A LONG-RANGE DIVISIBLE SANDPILE 


SUSANA FROMETA AND MILTON JARA 


Abstract. We study the scaling limit of a divisible sandpile model associated 
to a truncated a-stable random walk. We prove that the limiting distribution 
is related to an obstacle problem for a truncated fractional Laplacian. We also 
provide, as a fundamental tool, precise asymptotic expansions for the corre¬ 
sponding rescaled discrete Green’s functions. In particular, the convergence 
rate of these Green’s functions to its continuous counterpart is derived. 


1. Introduction 

In this paper we study the divisible sandpile model, which is a continuous version 
of the Abelian sandpile introduced in 1987 by Bak, Tang and Weisenfeld [1] as an 
example of a dynamical system displaying self-organized criticality. 

Up to our knowledge, the divisible sandpile was introduced by Levine and Peres 
in m- The main difference between the divisible sandpile and the Abelian sandpile 
is that instead of discrete particles of unit mass, each site can contain a continuous 
amount of mass. At the start there is an initial density of mass distributed on the 
lattice A lattice site is said to be a full site if it has mass at least 1. At each time 
step, each full site is toppled by keeping mass 1 for itself and distributing the excess 
mass among the lattice proportionally to the step distribution of a certain transition 
probability. As time goes to infinity the mass approaches a final distribution in 
which each site has mass less than or equal to 1. 

In [To] and Levine and Peres studied the scaling limit of this final distri¬ 
bution for the divisible sandpile in which the excess mass is split in every toppling 
equally among neighbors. In other words, they study the divisible sandpile model 
associated to the transition probability of a symmetric, simple random walk on 
Z"^. The proof is based within the framework of simple random walks. The corre¬ 
sponding scaling limit is related to the obstacle problem for the classical Laplacian 
operator. In m Lucas studied a divisible sandpile model that he calls the unfair 
divisible sandpile, in which the toppling procedure distributes mass to neighbors 
proportionally to a nearest neighbors, drifted random walk. The corresponding 
scaling limit of the final distribution is related to a so-called true heat ball. 

The so-called internal Diffusion Limited Aggregation model (iDLA) is a stochas¬ 
tic version of the Abelian sandpile introduced by Lawler, Bramson and Griffeath in 
[5| . In the mentioned paper, the authors studied the scaling limit of the random set 
of final occupied sites starting with m particles at the origin, which happens to be 
a ball with volume m in In m Levine and Peres also studied the scaling limit 
of the iDLA but starting from an arbitrary distribution of mass, they proved that 


2010 Mathematics Subject Classification. 60K35, 31C20, 35R35. 

Key words and phrases. Divisible sandpile. Green’s functions, obstacle problem, o;-stable laws. 
S.F. and M.J. were partially supported by CNPq and FAPERJ. 


1 



2 


SUSANA FROMETA AND MILTON JARA 



(a) 3400 iterations (b) 26500 iterations 


Figure 1. Truncated a-stable sandpile starting with mass 700 and 1300 
in two different sites. In each picture we use a scale of color that indicates 
the amount of mass at each site. 


considering the same initial distribution in the divisible sandpile and in the iDLA, 
the scaling limit of the final set of occupied sites in both models happens to be the 
same. In [12], a drifted version of the iDLA was also studied and the scaling limit 
was also proved to coincide with the scaling limit of the drifted divisible sandpile. 

Regarding the Abelian sandpile model some recent advances correspond to Peg- 
den and Smart, they proved in m that the Abelian sandpile converges in the 
weak-* topology to a limit characterized by an elliptic obstacle problem. Later 
in |9|, Levine, Pegden and Smart identified Apollonian structures in the scaling 
limit of the Abelian sandpile. Structures of the limit of the Abelian sandpile were 
previously observed in M- 

Our aim is to study a family of divisible sandpiles, which we call truncated, a- 
stable divisible sandpiles, that distribute mass not only to nearest neighbors. We 
present a scaling limit for the final distribution of a sequence of divisible sandpiles on 
which the excess mass is distributed proportionally to the transition probability of a 
truncated a-stable random walk. Figure [^represents a simulation of the truncated 
a-stable sandpile starting with mass in two different lattice sites. 

Green’s function estimates for the truncated a-stable random walk are a key 
tool in this work. In the context of simple random walks, such estimates are well 
known. A comprehensive reference, which also shows the power of these estimates 
as a tool in a variety of contexts, is Lawler’s book [^. In [^, Le Gall and Rosen 
obtained powerful estimates for Green’s functions of random walks in the domain 
of attraction of stable laws. 

In the specific context of truncated a-stable random walks with a € (1,2) in 
dimension d = 2, we present in Theorem |1.3| a precise asymptotic expansion for 
the corresponding rescaled discrete Green’s functions. In particular, the conver¬ 
gence rate of these Green’s functions to its continuous counterpart is derived. This 
asymptotic expansion is the most technical part of this work, and it is of inde¬ 
pendent interest. In particular, the optimal convergence rate is obtained only 
after the introduction of a non-trivial, nearest-neighbor correction to the transition 
probability of the truncated a-stable sandpile. Without this correction term, the 
convergence rate can be checked to be 
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Let /9 be a non-negative function on with compact support. We will use this 
function p to choose the initial mass density of the truncated a-stable sandpile. As 
the lattice spacing goes to zero, we run the sandpile with initial distribution given 
by p restricted to the lattice. Our main result, Theorem 1.2 states that the final 
distribution of the rescaled sandpile converges and that its limit is related to an 
obstacle problem for the truncated fractional Laplacian. Obstacle problems are a 
current subject of interest in the literature, see for example m, i. However, up 
to our knowledge the obstacle problem we need to deal with in this work has not 
been considered in the literature. 


1.1. Notations and results. The next few paragraphs are devoted to describe the 
random walk we use in the toppling procedure of the truncated a-stable divisible 
sandpile model. 

Let us introduce the following notations: 

For i? > 0, we denote Bn = {a; S |a;| < i?}. For a domain A in we denote 
A'-' = An 1?. We define the discrete boundary of the domain A, as 

dA'' = {x e 1?-, [x — (fl A and [x — ^,x + ^]'^ n A ^ 0}. 

We fix a parameter a S (1,2), and two positive parameters r < M. We define 
Ar,M '■= Bm \ Br- Now we will consider a function = Fi^r,M : —>■ R, which 
is essentially the indicator of the set A^ but for computational convenience, our 
function Fi takes specific values in dA]^j^. We define Fi{y) as the proportion of 
the square [y — ^,y + 5 ]^ that is inside the set that is 

Fi{y) = M([y- n a^^m), for y g z^, 

where p. represents the Lebesgue measure in R^. Notice that Fi is equal to 1 in 
^r,M \ equal to 0 outside A),)^ U OA^j^j. 

We define the probability jump pi : Z^ —)• [0,1] as 

Pi(y) = Ci(fcil{|y| = i} + foryGZ2\{0} (1.1) 

and pi(0) = 0, where ci is the normalizing constant and ki is a positive constant, 
which depends on a, r and M, defined as 




T—^dy — 




Fi{y) 

W 


( 1 . 2 ) 


We choose r and M in order to have ki > 0, see Lemma [4.1[ Then the probabil 


ity (1.1) is well defined. 


In order to obtain the desired convergence of the Green’s function, instead of 
working with a pure truncated a-stable random walk, we add a simple random walk 
multiplied by a carefully chosen constant. This is the constant fci defined in (1.2). 


This strategy leads us to cancellations that allows us to obtain convergence of order 
as Theorem 1.3 states. In particular, the convergence stated in Theorem 1.3 
does not hold without this correction. 

Now we run in lA our sandpile model starting with initial distribution p. When 
time goes to infinity, in Proposition |2.1| bellow we will prove that if every full site is 
toppled infinitely often, then the mass converges to a limiting distribution on which 
each site has mass taking values between 0 and 1. 
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Two individual topplings in the truncated a-stable sandpile do not commute, but 
the sandpile is Abelian in the sense that the limiting distribution does not depend 
on the ordering of topplings. This is the subject of Proposition |2.2| 

In the truncated a-stable divisible sandpile, as in the classical sandpile model, 
we use the odometer function as the fundamental tool to identify the limiting mass 
distribution. The odometer function measures the total amount of mass emitted 
from the point x G counted with repetitions: 

u{x) = total mass emitted from x. (1-3) 

Also as a consequence of the Abelian property, the odometer function does not 
depend on the order of the topplings. The total mass received by site x from other 
lattice points throughout the toppling procedure, is equal to udB'- Pi 

y), hence 

Lm,iu{x) = v{x) - p{x), (1.4) 

where p and v are respectively the initial and final distribution of mass, and Lm.i 
is the discrete Laplacian related to pi given by 

LMpf{x)= Pi(y)(/(ai + 2/) -/(ai))> (1-5) 

yeB-^jUdB-^ 

for a function / in Z^. 

We say that a function / is superharmonic at a point x G with respect to 
Lmb, if LM,ifix) < 0 . 

To study the limiting distribution we will use the following approach. First 
we construct a function 7 on Z^ whose discrete Laplacian is equal to p — 1, this 
will be done in Section Assuming we have this function, we consider its least 
superharmonic major ant 

s{x) = inf{/(x); Lm,i/ < 0 on Z^ and / > 7 }. ( 1 . 6 ) 

The majorant s is the solution of the discrete obstacle problem for the operator 
Am.i, with obstacle 7. 

The following lemma gives us an expression for the odometer. 


Lemma 1.1. Let p be a nonnegative bounded function on Z? with finite support. 
Then the odometer function for the truncated a-stable sandpile started with mass p 
satisfies 


u = s — "f 


where 7 satisfies Lm,i 7 = P — 
majorant of^. 


1 and s given by ( 1 . 6 ) is the least superharmonic 


Proof. The odometer function is nonnegative, thus u-l -7 > 7 . Since < 1, by (1.4) 
we have Lm,i(u -|- 7 ) < 0, hence s < u -|- 7 by (1.6). 

For the converse inequality let / > 7 be any superharmonic function lying above 7 , 
then Lm,i(/ — 7 — m) < 0 in the domain D = {x gZ"^] v{x) = 1} of fully occupied 
sites. Outside D, since u vanishes f — j — u > 0 and hence, using a maximum 
principle for Lm,i, we conclude that it is nonnegative everywhere. □ 


The expression obtained for u in Lemma |1.1| is called odometer’s equation. 
Lemma o gives us a way to find the odometer function once we have a func¬ 
tion 7 which satisfies Lm,i 7 = P ~ 1- This function is called the obstacle function. 
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Consider the random walk {5'm}m>o starting at 0 whose jump probability is 
given by pi, we define its Green’s function Gm,i as 

N 

GM.i(a;) = lim 

AT—)-oo 

m—0 

where = P[S'm = a^] and 

<^M,i = Lm,i\^\^ = pi(y)lpl^- (1-8) 

yeB-^UdB-^ 


Notice that i is the variance of the random walk with transition probability 
Pi- Sometimes in literature the function Gm,i defined in (1.71 is called potential 
kernel, see for instance [S] or [7]. Using the Markov property we can check that 
Lm,iGm,i{x) = 0 for a; 7 ^ 0 and Um,iGm,i(0) = —1. 

The obstacle function 7 can be, for instance, be equal to 


'y(x) = -- 


'MA 


E 


Gm,i{x — y)p{y) for all x 


(1.9) 


1.2. Scaling procedure. We will introduce a notation for the transition from the 
Euclidean space to the rescaled lattice 

For X S we write x"’” for the nearest point in -1? breaking ties to the right, 
that is x"’” = (x - ^,x + ■ 

For X S we write x"’'^ = ^ + 

For a function / in we write for its restriction to the lattice, that is 

/"- = /Uz- 

For a function / in ^Z^, we write for its extension to the Euclidean space, 
defined as /"’'^(x) = /(x”’”). 

For a domain A C we write = AH ^Z^. 

If A is a domain in -Z^, we write = A + 

For a domain A in K^, we define its discrete boundary with respect to the rescaled 
lattice -1? as 

n 

dA^’-'- = {x G iz2; X"’° (/ A and x"'° n A 7 ^ 0}. 


In the above notations, the index n will be omitted whenever is clear from the 
context that we are making reference to rescaled lattice ^Z^. 

Recall the notation \ Bn. We run our sandpile model in ^Z^ 

distributing the excess mass at each toppling, proportionally to the step distribution 
of the random walk with jump probability : ^Z^ —)■ [ 0 , 1 ], given by 


Pn{y) = Cn{kr 




Fn{v) 


’■^v\=A 


?2+“|y|2-|-a 


-) for y&-I?\ { 0 } ( 1 . 10 ) 

'/ n 


and Pra(O) = 0 , where c„ is the normalizing constant and the function Fn is defined 
as 

Fn{y) =n^ X p{y° n Ah^m) for y G -I?. (1.11) 

n 

Recall that p denotes the Lebesgue measure in . Notice that is equal to 1 
in At \ dAi and equal to 0 outside At U 9At 
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The constant c„ is the normalizing constant, and the constant is equal to 


,2-a 


kn — 


1 ,_^ Pniy) \ 

\y\a y |y|a )■ 


yGBjLU9Bii 


\y\^ 


( 1 . 12 ) 


As we mentioned before, we choose r < M to ensure that ki defined in (1.2) is 


positive. The existence of such constants is the subject of Lemma [4.1[ Moreover, 
as a consequence of the proof of Lemma |4.1| we have that kn is also positive for 
the same choice of the parameters r and M. Note that the truncation remains 
macroscopically constant. 

For a function / defined in define the discrete, truncated fractional Lapla- 
cian as 

LM,nf{x) = n°‘ Pn{y){f{x + y)- (1.13) 

Suppose we start the sandpile in with initial distribution given by the function 
Pn = Let us define the obstacle 7 „ as a function satisfying 


LM,n"fn{x) = p„{x) -1. 
The obstacle function can be defined for instance as 


Jn{x)=--^ - \ Y GM,n{^-y)Pn{y), 


(1.14) 


(1.15) 


'M,n 




where „ = LM,n\x\'^ and GM,n is a function which satisfies LM,nGM,n = —n^So- 
We call GM,n the discrete Green’s function. 

In Section we discuss how we define the Green’s function Gm.k and give an 
integral expression for it, that we anticipate here without proofs. All the details 
can be found in Section]^ For a fixed xq € with |a:o| = 1, define 


GM,n{x) = 


1 


cos(0 • x) — cos(0 ■ Xo) 


(271)2 


de, 


(1.16) 

(1.17) 


f [—n7r,n7r]2 i’M,niO) 

where 

■*/'M,n(6') = n“ Y Pniy)^^-cos{e-y)). 

yGB'^^UdB'^ 

The least superharmonic majorant for the obstacle 7 „ is a function s„ in 
defined as 

Sn{x) = inf {f{x); LM,nf < 0 on ^Z^ and / > 7„}. (1.18) 


By analogy with (1.3) we define the odometer function starting with mass pn in 
-Z2 as 

n 

Unix) := ^ X total amount of mass emitted from the site x. (1-19) 

It is easy to see that the odometer’s equation also holds, that is 

Un — 7n, ( 1 -^^) 


and the final distribution of mass Vn starting the sandpile with initial distribution 
of mass Pn is given by 

l^n ~ Pn T Gj^ n'^n- 
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1.3. Main results. The odometer’s equation (1.20) allows us to formulate our 
problem in a way that translates naturally to the continuum. Let c be defined as 
c := lim„_).oo Cm where c„ is the normalizing constant in (1.10), in Section|^we will 
prove that this limit exists and the convergence is of the order of n““. 

The continuous counterpart of the discrete, truncated fractional Laplacian Lm.k 
defined in (1.13) is the truncated fractional Laplacian whose integral form is given 
by 

f{x + y) + f{x-y)-2f{x) 


LMfix) = ^ / 


-dy, 


( 1 . 22 ) 


defined for a twice differentiable, bounded function /. 

In Proposition 4.3 bellow, it becomes clear why Lm is the natural extension of 
LM,n to the continuum. 

Let us fix a function p in with compact support, which represents the initial 
mass density. We dehne the continuous obstacle function as 


j{x) = - 2 “ “ / ^Mix - y)p{y)dy, 


(1.23) 




where = c jd^dy, and Gm is the Green’s function of the truncated frac¬ 
tional Laplacian defined as 


GMix) = 


cos{9 ■ x) — cos(d • xo) 


(2 t)^ 


where 


V’m(^') = c 


i’MiO) 

1 — cos (0 • y) 


de, 


dy 


(1.24) 


(1.25) 


and xq G is arbitrary but hxed. 

The idea behind the definition of the function Gm as the Green’s function of the 
truncated fractional Laplacian Lm, is that Gm is, in some sense, the inverse of the 


operator Lm- Such property will be discussed in Lemma 4.4 

We say that a function / on is superharmonic with respect to Lm if LmI < 0. 
Since we will often be working with functions that are not twice differentiable, it is 
convenient to define the truncated fractional Laplacian in a more general setting. 

We extend Lm by duality to a large class of distributions. A function / S Lj^^ 
is a linear functional in the dual of the space of infinitely differentiable functions 
with compact support G^. The symmetry of the operator Lm allows us to define, 

for / e 

{LMf,4>) = 

where (p G G^. This definition coincides with the previous one in the case f G G^. 

We say a function / G is superharmonic in an open set if for every 

nonnegative test function p G G^(R^) whose support is inside LI, we have 

{f,LMp)<0- (1.26) 


Throughout the text, superharmonic function will always mean super harmonic 
with respect to Lm- 

Now we consider the least superharmonic majorant 

s(x) = inf{/(x);/ is continuous, superharmonic and / > 7 }. (1-27) 

The majorant s is the solution of the obstacle problem for the truncated fractional 
Laplacian with obstacle 7 . 
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The odometer function for p is given by 

u = s — (1-28) 

We say that a sequence of functions /„ S converges to a distribution T in 
the weak-* topology if for any test function </> S we have 

lim = {T,4>). 

n—^oo 

Our main result states the convergence of the final distribution of the truncated 
a-stable sandpile with respect to the weak-* topology. 

We denote by the space of the twice differentiable functions with compact 

support. 


Theorem 1.2. Let a € (1,2). Consider p G as the initial density of the 

truncated a-stable sandpile. Let Vn the final distribution of the sandpile in 
starting with initial distribution pn = p". Then 

P + Lmu 


in the weak-* topology, where u is the limiting odometer function defined in (1.281. 


As we mentioned before, we obtained the next following convergence result for 
the Green’s functions. 


Theorem 1.3. Let Gm.u o.'rid Gm be respectively the rescaled discrete Green’s func¬ 
tions and the continuous Green’s function defined in (1.16) and in (1.24). Then 


there exists a constant K > 0 which only depends on a, such that for all x 0 and 
all n, holds 

Gm (a;) - Gm, n (a:)+/3„ log I a; I < ^(l+ (1.29) 

where fin is an explicitly computable constant (given in terms of the probability 
jumps) which converges to zero at rate n““. 

In particular Gm.u converges towards Gm uniformly on compacts o/K^ \ {0}. 

1.4. Organization of the paper. In this work we will follow the next structure 

convergence of Green’s functions ^ convergence of obstacles 

=> convergence of majorants 
convergence of odometers 
convergence of final distribution. 

In Section|6]^ treat the convergence of the Green’s functions, which is the subject 
of Theorem Ol In Sectionj^we prove the convergence of obstacles as a consequence 
of Theorem 1.3| Although we do not use all the information that Theorem [T^ gives 
in the proof of the convergence of obstacles, we do need such a strong result for 
the convergence of majorants. In Section]^ we prove the convergence of majorants 
and the convergence of the odometer functions. A fundamental tool for Section 
1^ is the continuity of the majorant ( 1.27[ ) which is the subject of Section]^ In 
Section we obtain our main result, Theorem |I.2| In Section we put together 
some important analytical estimates that we will use throughout this work. Proof 
of various auxiliary results are in Sections and El 
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2. Existence of the odometer function and Abelian property 


We start with an initial distribution of mass p in 1? with finite support. At each 
time step we topple a full site, recall that a toppling in the site x consists in leaving 
mass 1 to cc and distribute the excess mass among the lattice proportionally to the 


step distribution of the random walk (1.1). 


In this section we prove that if every full site is toppled infinitely often, the mass 
distribution converges to a limiting distribution on which each site has mass less 
than or equal to 1. We also prove the Abelian property which states that the final 
distribution does not depend on the order of the topplings. 

We define a toppling scheme T as an infinite sequence of indexes in 1? in which 
each full site that is initially full or becomes full after the realization of the previous 
toppling appears in the sequence infinitely often. We also say that a toppling scheme 
is legal if it only topple full sites. 

Let uj be the mass distribution after the toppling of the first k sites listed in T, 
and (x) the mass emitted from x up to the fc-th toppling. 

Proposition 2.1. Suppose we start with an initial configuration p with finite total 
amount of mass and bounded support. Consider a legal toppling scheme T for this 
initial configuration. Then as k goes to infinity, and tend to limits v and u. 
Moreover the limiting configuration v satisfies v < 1 inlP' and u = p + Lm,iu, for 
Lm,i defined in (1.5). 


Proof. We claim that there exists a bounded set K which contains all the possible 
sites with positive mass throughout the realization of the topplings. To see that 
denote by S the support of the initial distribution p and let m be the total amount 
of mass. If a lattice site y outside S becomes full at certain instant of the toppling 
scheme T, then there exists a sequence yi,y 2 , - ■ ■ ,yt. of full sites such that yi = y, 
yi G S and j/i+i G {yi\ + Note that, since the mass is preserved, £ < m. 

Then all possible full sites are inside the set K := S + i?m(M+i) • This allows us 
to conclude that the final set of sites with positive mass is contained in the set 
K := K + Bm+1- 

Let us consider the function 


where ^ was defined in (1.8). 

Note that Wk is uniformly bounded. 

Assume that x is the /c-th site to be toppled, then the mass at point x is modihed 
by j3k{x) = Vk-i{x) — Ufc(x), and this amount will be transfered to the lattice 
according to the toppling rule. Hence, 


Wk{x) - Wk-iix) = 


{My) - Vk-i{y))'^^ 




^M,l 


= -/3fc(ai)- 


'M.l 


Y ^k{x)pi{y - x) 

y^x 




'MA 


= /3k{x)LM,i- 

( 

= 


^ MA 
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Since uu is the sum up to k of all the relevant we get: 

Wk = Wo+J2 ^kix). 

Then for every x, Uk{x) is a bounded increasing sequence, then it converges to 
certain a function u{x). 

The relation Vk{x) = p{x) + LM,iUk{x) holds for all finite time fc, so the conver¬ 
gence of Vk is a consequence of the convergence of Uk- Moreover, its limit v satisfies 
V = p + Lm,iu. 

Finally a point x is either never toppled, in which case we have Vk{x) < 1 for 
all k, or it is toppled infinitely often, and then when a toppling occurs at time k at 
the point x we have i^ki^) < 1. In both cases we conclude that the limit v satisfies 
v{x) <1. □ 

The next proposition states the Abelian property for the truncated a-stable 
sandpile. 

Proposition 2.2. Consider two legal toppling schemes Ti and T 2 for an initial 
configuration p. For i = 1,2, let and be respectively the final distribution 
and the odometer function for the toppling scheme Ti. 

Then and . 


Proof. Assume that at time k we topple the site x^ in the toppling scheme Ti. We 
will prove by induction that u'^^{xk) > u^^(xk). 

This property is trivially true for A; = 1. Suppose that this holds for all i < k. For 
X Xk we have two possibilities: either x is not toppled before time k in the scheme 
Ti, in which case {x) = 0, or x is toppled before time k. In this case we consider 
the last index i in which the site x is toppled. Then > uj^{x) = u^^{x). In 

both cases 

u^^{x) >ul^{x). ( 2 . 1 ) 

Since Ti is a legal toppling scheme, 

< I < ( 2 . 2 ) 


The initial configuration is the same for both toppling schemes, so equation (2.2) 
says: 

LM,lU^^{Xk) < LM,iul^{xk). 

That is, 

U^^{xk) - U^^{xk) > y] Pl{x - Xk)(u^‘^{x) - Mfc^(x)). 


Each term in the sum above is nonnegative by (2.1), then we conclude that 
li^^fxk) > M^^(xfc). Since the site Xk is toppled infinitely many times, it follows 
that . The same argument shows the reversed inequality, which concludes 

the proof. □ 


3. The obstacle function 


Equation (1.20) gives us a way to find the odometer function for the truncated a- 
stable sandpile, provided that we have a function 7 „ that satisfies condition (I.I4). 
This section is devoted to the construction of our obstacle function 7 „. 
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Recall the definition of Pn given in (1.10 1 and define the jump probability on 
as: 

Pniy)=Pn{^)- (3.1) 

For a function / on Z^, we define the generator Fm.ti of this random walk, as: 

LM,nf{x)= Pn{y){f{x + y) - fix)). (3.2) 

For each n we consider in If’ the random walk {Slf^}m>o starting at the origin, with 
law given by p„. Define p^ix) = PI'S”))! = x\. The function Gm,™ given by 


N 

Gm, nix) = ^lil^ Y [Pn (0) - Pnix) 

m —0 

verifies LM.npM. njx) = 0 for cc ^ 0 and TM,nGM,n(0) = — 1 . 
The function ( |3.3[ ) can be written as (see my- 

f cos (0 ■ x) — I 


(3.3) 




(27r) IpM.niS) 


d0, 


(3.4) 


where ipM.ni^) is one minus the characteristic function of the probability of (3.1), 
that is 

'0M.n(6') = Y Pn(2/)(1 - cos( 6 < • ?/)). (3.5) 

We are not actually interested in Gm.u, but in the value of its Laplacian, then we 
can add constants in a convenient way in order to obtain convergence when the 
function is properly rescaled. So we fix some xq in such that |a:o| = 1 and notice 
that the function 


GM,nix) = 


1 


(271)2 


cos (0 • x) — cos (0 • nxo) 

i’M,ni0) 


de 


(3.6) 


also satisfies LM,nGM,nix) = 0 for a: 7 ^ 0 and TM,nGM,n( 0 ) = — 1 . 

We define the function Gm.u '■ ^Z^ —>• M given by 

Gm, nix) = rf~°'GM,ninx) (3.7) 

This function satisfies LM,nGM,nix) = 0 for a; 7 ^ 0 and LM.nGM.niO) = —'rf. A 
change of variables shows that Gjvf.n has the integral form given in (1.16). We call 
GmjTx the Green’s function associated to the operator LM,n- 
Define 

GM,n*Pnix) = ^ Y. GM,nix - y)pniy)- 

The operator Tm.m commutes with the sum, therefore 

LM.niGM.n * Pn)ix) — iLM,nGM,n) * Pnix) = Pnix). 

Now we are ready to define the obstacle function 7 „ on as 


Inix) = - 2 - Gm.u * Pnix), 


(3.8) 


' M.i 


where cr^n = LM,n\x\^, as defined in (1.15). 
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4. Analytical estimates 


In this section we put together some of the analytical estimates that we use 
throughout this paper. We also state some qualitative properties of the Green’s 
function, the obstacle function, and the odometer function. The proofs of these 


propositions can be found in Section 10 


Recall from (1.2) the definition of ki. The next lemma states that we can choose 
the parameters in the definition of ki in order to ensure that ki is positive and 
consequently the probability jump pi is well defined. 

Lemma 4.1. There exist constants r < M such that ki is positive. 

The next lemma tells us a simple fact concerning the normalizing constants c„ 
of the jump probability Pn (recall from (l.lOl), and the constant c = lim. 


n—>-oo Gn • 


Lemma 4.2. There exists a positive constant C which depends on a, r and M 
such that 

C 

|cn — cl < — for all n. 

In the next lemma we show that for a sufficiently smooth function / defined 
in IR^, its discrete fractional Laplacian Lm.u defined in (1.13) approximates its 
fractional truncated Laplacian Lm defined in (1.22). 


Lemma 4.3. If f : 
set. Then 


^ K IS twice differentiable in and K CM.'^ is a compact 
\LMf{x) - LM,nf{.x)\ < 


for X € K, where the constant C depends on a and on the first and second partial 
derivatives of the function f in K + Bm • 

As we mentioned before, the Green’s function Gm (recall the definition from 
(1.24)), is in some sense the inverse of the truncated fractional Laplacian Lm- This 
is the subject of the next lemma. 

Lemma 4.4. Let </) : — 

support. Then 


a infinitely differentiable function with compact 


Lm{Gm *</>) = —(f- 

The next two propositions give us a better understanding of the behavior of the 
Green’s function Gm defined in (1.24). 


A change of variables shows that 


GM(x) = ^(Gi(^)-Gi(t)). 


(4.1) 


We have the following estimates for the function Gi 


Proposition 4.5. There exists a bounded function g in Bi, and a bounded function 
h in Bf, so that 

„,2 1 


Glix) = 


:(27r)2 


+ g{x), for all |a;| < 1. 


and 


Gi(x) = -^ log |a;| + hi(a;) + 


X 


2-a 


for all |x| > 1, 


(4.2) 


(4.3) 
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where hi is the bounded function in Bf, given by hi{x) = dQ, 

where uj is any unitary vector. 


We have the following estimates for Gm- 


Proposition 4.6. There exists a positive constant C independent of x and M, 
such that 


Gm{x) 


1 

c(27r)2 |xp““ 


c(27r)2 


C 

M 2 -«’ 


for all X S Bm 


(4.4) 


and 


Gnix) 


TTCr 


■ log \x\ - hM{x) - Sm 


M 


< 


G 


— for all X G B 

— rv *' 


M- 


(4.5) 


where 6 m is an explicitly computable constant, converging to o,s M goes 

to infinity, and hM is a bounded function in Bfj, given by the formula hM{x) = 
Ib i^i \Bi |e|2 ^ d0, where uj is any unitary vector. 

“ -ff 

The function G(x) := Green’s function of the fractional Lapla- 

cian whose integral form is given by (see m) 


Lf{x) 


f{x + y) + f{x-y)-2f{x) 

\y\2+a 


dy, 


(4.6) 


for a twice differentiable, bounded function /. 

Also note that the function-Vloglxj is the Green’s function of a Brown- 

ian motion (see 0). then Proposition |4.6| states that the Green’s function Gm is 
an interpolation between the Green’s functions of a Levy a-stable process and a 
Brownian motion. 

The corollary bellow is a direct consequence of (4.4). 


Corollary 4.7. The function Gm{x) converges to G{x) — uniformly on com¬ 

pacts sets o/K^. 


In some situations we will need to integrate the first partial derivatives of the 
Green’s function Gm over finite sets. In the following proposition we justify such 
operations. 


Proposition 4.8. The first partial derivatives of the Green’s function Gm o-re 
locally integrable. 


Later on we will use that the odometer function u defined in (1.28) has compact 
support. The next lemma will be used to prove this fact. 


Lemma 4.9. There exists a ball 12 such that the continuous obstacle function 7 
defined in (1.23) is concave outside 12. 


Proposition 4.10. 

pact support. 


The continuous odometer function u defined in (1.28) has com- 
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5. Continuity of the majorant 


In this section we prove the continuity of the majorant s defined in (1.27). 


The continuity of the majorant, in the previous related works (see for exam¬ 
ple HD) is a simple consequence of the mean value inequality for the classical 
Laplacian operator. We can translate the same technique to our context and for 
that, in this section we present a mean value inequality for the truncated fractional 
Laplacian. We follow an idea from Caffarelli and Silvestre m. where they proved 
similar results for the fractional Laplacian. However, in our case, an extra care 
must be taken since we are dealing with a truncation parameter that may change 
when we do simple operations like change of variables on integration. 


Recall from (1.22) and (1.24), the definition of the truncated fractional Laplacian 
Lm and its respective Green’s function Gm- For A > 0, consider the Green’s 
function G m . 


We will construct a family of approximation of the identity. For every A > 0 
we stick a paraboloid from below to cut out the singularity of the function Gm at 
x = 0. In this way we obtain a twice differentiable function that we call F;^(x). The 
parameters of the paraboloid can be chosen in a way such that the function rA(x) 
coincides with Gm{x) when x is outside a ball of radius 1, and inside the ball Bi 
the function F;),(x) coincides with the paraboloid. 

Next we rescale F^, in the following way. 


rA(ai) = y^rA(f). 


(5.1) 


A change of variables shows that 

Gm{x) = A2-“Gm(Ax) - A2-“Gm(Axo), (5.2) 

so the function T;, coincides with Gm{x) — Gm(Axq) when x is outside Bx. 

Let us consider then the family of functions 

fA = rA + GM(Axo) (5.3) 


and note that Fa coincides with Gm outside Bx- Moreover Fa^ < rA 2 if Ai > A 2 . 

Proposition 5.1. For all A > 0, —Lm^x is a positive continuous function in L^. 
Moreover 


—LM^x{x)dx = 1 . 


Proof. We have the following rescaling property 

LMTxix) = LMfxix) = ^LMTxi^). (5.4) 

A^ ^ A 

Then a change of variables shows that it is suffices to prove that for all A > 0, 
—Lm^x is a positive continuous function in and — Lxix)dx = 1. 

We prove first the —LmTx is positive. 






SCALING LIMIT FOR A LONG-RANGE DIVISIBLE SANDPILE 


15 


If x ^ -Bi, then rA(x) = G m{x) and for every other y, rA(y) < GM{y), then: 

2rA(a;) - rA(x - y) - rA(x -h y) 


-LMT.ix) = 


\y\ 


2+q: 


-dy 


c 

> 

- 2 


= 0 


2Gm (x) — Gm [x — y) — Gm (x + y) 


\y\ 


2+q: 


dy 


since Gm is the Green’s function of Lm. 

A A 


If a; € Bi \ {0}, there exist a, xi ^ x and a positive 6 such that Gm (• — xi)-|-(5 
touches Fa from above at the point x. Then 

2rA(x) - rA(a; - y) - rA(a; + y) 


-LMT.ix) = 


\y\ 


2+a 


-dy 


c 

> - 
- 2 


= 0 


2Gm (x — xi) — Gm [x — y — xi) — Gm (x -\- y — xi) 

A A -A- 


\y\ 


2+q: 


since Gm is the Green’s function of Lm . In the second line above we cancelled the 

A A 

parameter 6. 


If a: = 0, then Fa attains its maximum at x: 

2FA(a;) - Fa(x - y) - Fa(x + y) 


-L„rx(x) = 5 / 


\y\ 


2+q 


dy > 0, 


because we are integrating a positive function. 

To show that — L\ix)dx = 1 we consider a smooth function w such that 
Iaj(x) < 1 for every a: € uj(x) = I for every x € Bi and suppw C B 2 . Let 
WR(a;) =a;(|), 

/ —LMFA(a;)dx — 1 = lim ( — Lm (Fa — Gm ), wb) 

Jr2 a ' ^ R^O ^ A ^ A ^ / 

= lim (Fa — G M , —L M wb) 

R^O ' A ’ X 

= 0 

since Lmlur goes to zero uniformly on compacts sets, and Gm —Fa is an function 
with compact support. □ 

Now let us define r/x = 

Proposition 5.2. The family yx is an approximation of the identity as A —> 0 m 
the sense that, if f is a continuous function in then 

Vx * f ^ f when A —> 0 uniformly on compacts. 

Proof. We will prove first that the collection of functions {LmFa,A > 0} is uni¬ 
formly integrable for A sufficiently small. That is, for all e > 0 there exists R > 0 
such that 




\LMrx{y)\dy < 
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for A sufficiently small. 

Fix y outside 82- Hence r>( 2 /) = GM{y), we have 


-LMTxiy) 


c f 2GM{y)-Txiy - z)-Tx{y + z) 

A / -^- dz. 

2 Jb^ 

A 


Since for y 7 ^ 0, LmGm (y) = 0, the equation above can be written as 


-LMTxiy) 


c r GM{y-z)- Txiy 

h + 12- 


z) , c 
-d, + 5 


Gm (y + z)- rA(y + z) 


\2+a 


dz 


Since the support of Gm — r>, is contained in Bi, the integral Ii above can be taken 
over the set Bm. H B{y, 1 ), and the integral I 2 can be taken in Bm n B{—y, 1). 

Since Gm is bigger than Fa in Bi, for a positive constant C independent of y, 
we write 


h 


< 


< 


cj 

JB(-y,l) 


G M (y z') 

|z|2+“ 


dz 


G 

|y|2+a 



Gm {z)dz, 


where the const ant G have changed from one line to the other. 
By Corollary 4.71 Gm converges uniformly on Hi to G — 
tioii of ■ 


c(27r)^ 


(where G is the 


Green’s function of the fractional Laplacian) when A goes to 0. 
Then we can pass the limit inside the integral: 


lim 

A->0 


Bi 


GM{y)dy=f {G{y) 


c{2tt) 


2 )dy- 


Since G is integrable in Hi, for A sufficiently small, there exists a constant G > 0 
independent of A and y such that 


h < 


G 

|y|2+a 


for y G B^. 


(5.5) 


The exactly same argument proves that (5.51 also holds for / 2 . Then we have 


-LMTxiy) < for all y G H^. 

Take R > 2, then 

f —L mL x{v)dy < f — 7771 — dy 

Jm'^XBr ^ Jts.’^xbr lyP "^" 

27rG 

which proves the uniform integrability. 

Now fix a compact set K, we will prove that t / as A goes to zero, 

uniformly on K, assuming that the function / is continuous. 
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Recall that ri\{x) = —Lm^x{x) = —^LmTx{^). We have, for a; G RT 
v\*f{x)-f{x) = [ v\{x-y)ifiy)- f{x))dy 

= f -LMTx{y){f{x - Xy) - f{x))dy 
= f -LMTx{y){f{x - Xy) - f{x))dy, 


(5.6) 


where in (5.6) we used that the support of LmTx is contained in this can 

be easily seen from the fact that Fa coincides with Gm outside B^. 

By the uniform integrability of {LmTx, A > 0}, we fix e > 0 and take i? > 0 so 
that 

/ —LMTx{y)dy < e, for all A sufficiently small. (5.7) 


Then we write, for x € K 
V\ * f{x) - f{x) = 


/ Bj^nB 


-LMrx{y){f{x - Xy) - f{x))dy 


(^+ 1 ) 


+ 


-LMrx{y){f{x - Xy) - J{x))dy 




—: Ji + J2 

Estimate of Ji: Since / is uniformly continuous on compact sets, we choose A 
small enough, so that 


Then 


\f{x - Xy) - f(x)\ < e, for all x G K,y G Br. 


|-^i|<e/ -LMrx{y)dy = €. 

JR2 ^ 


Estimate of J 2 ' For y G we have Xy G Bm+\ C Bm+i- Then in the 

integral J 2 , we are only considering / taking values in the set K + B(^m+i}- Since 
/ is continuous, in particular is bounded in RT + let say by A. 

Then, by (5.7), we have 

IJ 2 I < 2A / —LMTx{y)dy < 2Ae. 

Since e is arbitrary, we conclude our proof. □ 

In the next proposition we prove a mean value inequality for superharmonic 
functions with respect to the operator Lm- Recall from (1.26) the definition of 
superharmonicity for functions that are not twice differentiable. 

Proposition 5.3. Let f be a eontinuous function. Then f is superharmonic in an 
open set U if and only if 

f{x) > 7?a * /(x) 

for any x in U and any X satisfying B{x, A) C U. 


(5.8) 
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Proof. Suppose LmI < 0 in C7. 

Let r > 0 such that, r > Ai > A 2 , then — Laj is a nonnegative smooth 
function supported in Br- If Lm/ < 0 in B{x,r) then: 

Using the self-adjointness of Lm- 

(/, (-LM)rA,(- - x) - {-LM)fxA- - ^)) < 0. 

Therefore 

< (/,?7A2(--x)) 

< f*r]x^{x). (5.9) 

By Proposition |5.2[ since / is continuous, we have uniformly on compact 

sets, as A —0. 

We take A 2 —0 in (5.9), and this concludes our proof. 

The conversely part follows easily. □ 


In Proposition |5.3| we obtained the desired mean value inequality for the trun¬ 
cated fractional Laplacian. Before proving the continuity of the majorant s defined 
in (1.27) we need a few lemmas. 

Given a bounded open set U, define the function 

(x) = inf{/(ai) : / is continuous, superharmonic on U, and / > 7 }. (5.10) 

By Proposition |4.10[ there exists a ball 12 which contains the support of the odome¬ 
ter function u. We fix this ball 12 and assume that U satisfies 


n + B 


M+l 


c u. 


(5.11) 


Lemma 5.4. Let U be a bounded open set satisfying (5.11). Then the function s^ 
defined in (5.10) satisfies 

s^(x) >rix* s^(x), 
for all X € and A sufficiently small. 


Proof. Note that since the infimum in (5.10) is taken over a larger set than the 
infimum in the de finitio n of s, we have that s^ < s. 

By Proposition 4.10 s coincides with 7 outside the ball 12. Then 7 < < s, 

and hence, s^ also coincides with 7 in 12 °. 

Consider a continuous function / bigger or eq ual than 7 , and superharmonic in 
U. Then f > s^ and moreover by Proposition 5.3 for x S 17 and A sufficiently 
small we have 

/(x) > 77 a * /(x) > 77 a * s^(x), 
taking infimum over /, we have: 

s^(x) > 77 a * s*^(x) for X G U. 

By our choice of 17, the function 7 is superharmonic in 17° 


(5.12) 

Then, using Proposition 


5.3 for cc S 17° and A sufficiently small 

s^(x) = 7 (x) > 77 a * 7 (x) 


holds. 
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Notice that supp(77a) C -Bm-i-i for A < 1. Since U satisfies (5.11), we can replace 
7 by in the equation above, and obtain 

s^ix) > 7 a * s^ix) for X S t/°. (5.13) 

□ 


By (5.12) and (5.13) we conclude our proof. 


The next lemma states the continuity of . 


Lemma 5.5. Let U be a hounded open set satisfying (5.11). Then the function s 


defined in (5.10) is continuous. 


Proof. We have noticed before that coincides with 7 in 11°, then it is continuous 
in 11°. We will prove that it is continuous in U. 

By Lemma |5.4[ we have 

> 7a * s^, 

for A sufficiently small. 

Consider now a bounded, open set Ui, such that, Ui D U + Bm. 

Since 7 is continuous, by Proposition |5.2| we have that 7 a * 7 7 uniformly in 

Bi as A —?> 0, then for e > 0 we take A small, so that 


7 a * s^{x) > 7 a * 'y{x) > j{x) — e, for x G Ui. (5.14) 

Note that 7 a * is a continuous function. By Lemma 5.4 and Proposition 5.3 
we conclude that it is also superharmonic. 


Consider the function g = ma.x{r]\ * s^ + £, 7 ). By (5.14), g coincides with 
7 a * + c in Bi, and hence it is superharmonic in U. Since g is above 7 , we have 

g > s^. We have proved 

7 a * s^ix) < s^{x) < 7 a * s^{x) + e, for x G U. 

Thus s^ is an increasing limit of continuous functions in U, and hence lower semi- 
continuous. Since s^ is defined as an infimum of continuous functions, it is also 
upper semicontinuous. □ 


By Lemmas 5.4 and 5.5 and Proposition 5.3 we have that s^ is continuous and 
superharmonic. Since s^ is above the obstacle 7 , we obtain s = s^. 

We have just proved the following proposition which states the continuity of s. 


Proposition 5.6. The majorant s defined in (1.27) is a continuous function. More¬ 
over, if U is a bounded open set satisfying (5.11), then s = s^. 


6 . Convergence of Green’s functions 


In this section we prove Theorem |1.3| which tells us about the convergence of 
the Green’s functions. We refer to Section im for the proofs of the auxiliary results 
presented here. 


Recall from (1.16) and (1.24) the definition of Gjvi.n and Gm- The first step in 
order to prove Theorem 1.3 is to investigate the convergence of ipM,n to i/>m defined 
in (1.17) and (1.25) respectively. 


A simple change of variables shows that 

1 




( 6 . 1 ) 
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For 6 G i3j_, by (10.7) we have that 




for a bounded function qm ■ Bj_ 


In the same way, for 9 G , by (10.9) 


( 6 . 2 ) 


(6.3) 


= Cal^'r +hM{0), 

for a bounded function Hm ■ 

M 

This suggests that we must split the integrals that define Gm and GM,n in three 
different regions: [—rtt, nir]^ \ B_^ and \ [—nir,nir]^. Recall from (1.11) 

the definition of Fn. 

For 9 G [-n7r,n7r]2, denote ^mA^) = ■ Then 

we have 


'^mA^) = n‘^Cnkn X! + Cu^mA^)- 

\v\=k 

Note that for 9 G [—mr, nir]^, we have ^ is bounded, then we use Taylor’s Theorem 
to ensure that there exists a positive constant G, such that 


2(1 — cos( —)) + 2(1 — cos( —)) — 


92,, w ^C\9\ 


^ 4 ’ 


then 


i’MAA - Cn'0M.n(0) “ 


cAnW / C\9\ 


T 2 —a 


< 


(6.4) 


The next two lemmas describe the convergence of i/iM.n to i/’m when n goes to 
infinity. 


Lemma 6.1. For 9 G Bjl we have 


\'>PmA^) — -iPm{0)\ < — U -, 

c 


(6.5) 


where the constant G only depends on a. 

Lemma 6.2. For 9 G [—n7r,n7r]^ \ B±_ there exists an uniformly bounded family 
of functions such that 


V:m.„(0) = -V’m( 0) + 
c 


( 6 . 6 ) 


Moreover there exists a positive constant C depending only on a such that sequence 

r„ satisfies 

(i) ^A{9)\9A’^) < C|0|i+“ for 1 = 1,2 

(ii) A(r„(0)|0|2+“)<C|0|“. 

Let us introduce the following notations: 

cos(0 • x) — 1 


GM^niA) — 

Gm{x) = 


Ibi i’MAA 

'M 

f cos{9 • t) — 1 


IB 


flMiO) 


d9 + 


d9 + 


f cosifi ■ x) 

[—n7r,n7r]^\B i 
IfT 

f cos(0-a:) ^^^ 
iiM{S) 


d9 (6.7) 


( 6 . 8 ) 
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The next lemma states the convergence of the auxiliary functions defined in (6.7) 
and (6.8). In the sequence we will conclude the proof of Theorem 1.3 which is the 
subject of this section. 


Lemma 6.3. There exists a eonstant C which only depends on a, such that 

1 1 
I 'T*l I 'T- I 2 , 


c ~ ~ C 

\—Gm{x) - GM,n{x)\ < — 


1 


Proof of Theorem |1. Note that 
GM,n{x) = 


1 


;{GM,n{x) — Gm, n{xo)) 


(271)2' 

Gm{x) = ^(Gm(x)-Gm(^o))- 


By Lemma |6. 3 


(27r) 


— Gm{x) - GM,n{x) 


G 1 

— — (^ 1 —r 

71^ \ \x\ 


tp)- 


(6.9) 


Let us define Pn = — l), and note that /3„ converges to zero with rate 

due to Lemma 


Theorem 11.3 


4.2 


Finally, using Proposition 


4.6 


we conclude the proof of 

□ 


7. Convergence of obstacles 
Fix X G and y G define 

0-n{y)-=\GM,n{x"-y)-l GM{x-z)dz. (7.1) 

jyU 

Lemma 7.1. The function orf satisfies 

Kfa)l s 

where f is a positive function in and G is a positive constant which depends 
on a. 

Proof. Write 

(^n{y) = \{GM,n{x'-y)-GM{x''--y))+ [ {Gm{.x-''- y)-GM{x-z))dz. (7.2) 

JyO 

For z G 7/'^, by Taylor’s Theorem, there exists a positive constant G such that 

\Gm{x''' -y)- Gm{x - z)\< —\VGm{x'''' - y)\, 

n 


Integrating we obtain 


/ \Gm{x''’'- y) - Gm{x - z)\dz < 

Jv^ 


G|VGm(x”-2/)| 


(7.3) 


On the other hand by Theorem |1.3| there exist a constant G > 0 depending only 
on a, such that 


By Proposition 4.8 the function |VGM(a;)| is in 


\Gm{x'-'- -y)- GM,n{x-- - y)\ < ^[/(x''' - y) + |— 


(7.4) 
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where / is a function in 

The function is not locally integrable. To overcome this difficulty, since x'' 
and y are both in -1?, we write 

^ n ’ 


C-- - y\ 


2 < 


C-- -y\' 


This last inequality, together with (7.2), (7.3) and (7.4) conclndes the proof of the 
lemma. □ 


Proposition 7.2. Let "fn andbe the obstacles defined in (3.8) and (1.23). Then 
as n goes to infinity 

7 ^ —>■ 7 uniformly on compact sets o/K^. 

Proof. We have „ —>■ We only need to prove that {GM,n * Pn)^ —t Gm * P 
uniformly on compact sets of 
Since pn = p " we write 

(Gm.u * p)^ix) — Gm * pix) = Gm.u* p{x") — Gm * p{x) 

= p{y)GMAx''' -y)-I p{y)GM{x-y)dy 

= 51 PWl<kv), 


where affy) is given by (7.1). 

Recall that p has compact support, let us say E. By Lemma [7T| there exists a 
positive locally integrable function / which satisfies 

C 


\GM,n* p{x'') - Gm * p{x)\ < 


il+ct 


51 -y) 


yeE--- 


< 


G 


f{x - y)dy, 


IE 


where the constant G has changed from one line to the other. Since / S and 
a > 1 , we conclude the proof. □ 


8 . Convergence of odometers 


In this section we prove the convergence of majorants, and we obtain as a con¬ 
sequence the convergence of odometers. Before that we need a few lemmas. 

Recall from (1.15), (1.18) and (1.19) the definition of the discrete obstacle func¬ 
tion 7 „, the discrete majorant s„ and the discrete odometer function and recall 
the definition of its respective continuous counterparts 7 , s and u from (1.23), (1.27) 
and (1.28) respectively. 


Lemma 8.1. There exists a ball LI which contains the support of the odometer 
functions. That is 

supp(u) C LI and (^supp(tt„) C LI. 

n 
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Proof. By Lemma |4.9[ there exists a ball such that the obstacle function 7 is 
concave outside ft. In Proposition |4. 10| we proved that this ball ft also contains the 
support of the odometer function u. 

We can write the discrete obstacle function 7 „ as 

IxF 

'yn (^) — 2 {G M.n M) * (^) (^) ■ 

^M,n 

We use Theorem |1.3| and repeat the argument used in the proof of Lemma |4.9| to 
conclude that there exists a ball (bigger if necessary) such that the obstacles 
functions 7 „ are concave outside ft. 

Then, the concave envelope r„ of the obstacle 7 „, defined as the least concave 
function lying above 7 „, satisfies that r„(a:) = ')n{x) for all x outside ft and for 
all n. The same argument used in the proof of Proposition |4.10| concludes our 
proof. □ 


Lemma 8.2. Let 7 „ and Sr 


be the obstacle function and the majorant given by 
is a bounded set like the one in Lemma 


( |1.15| ) and ( |1.18| ). If ft C 

Sn{x) = inf {/(cc); Lm,™/ < 0 on ft' and f > 7 „} 


11 


then 


( 8 . 1 ) 


Proof. We observe first that LM,nSn = 0 on the set = 1}, since 

I^M.n^n — I'M.nip^n T ^n) — Pn “t” 1 — Pn 1- 

Let US call 

s)J(a:) = inf {f{x) : LM,nf < 0 on ft' and / > 7 „}. 

We want to prove that = s„. Since the infimum is taken over a strictly larger 
set, the inequality s^ < Sn is trivial. 

For the converse inequality consider a function f > jn such that LM,nf < 0 on 
ft '. Since Lm.uSu = 0 on D„, the difference f — Sn satisfies LM,n{f — Sn) < 0 
on Dm then by the maximum principle f — Sn attains its minimum at a point x 
outside Dm where Sn{x) = ^n{x). Since f > Jn conclude f > Sn on ft ' and 
hence everywhere. □ 


Recall that Lemma l4~4l states 

Lm{Gm * 4>) = —4>, for a function (j) e Gf°. (8.2) 

The next lemma studies the equation above removing the smoothness assumption. 

Lemma 8.3. Let f be a bounded, measurable function, with compact support. Then 
Lm{Gm * /) = —f 'In the sense of distributions. 

Proof. For a function (j) G Gf°, holds Lm{Gm *</') = Gm * {Lm4>)- We have 
{Lm{Gm * f),(t>) = {f,L m{Gm * f)) 

= (/,-</>) 

= i-fA), 


where in the second line we used equation ( 8 . 2 ). 


□ 


Recall from (1.26) the definition of superharmonicity for functions that are not 
twice differentiable. The next lemma is a consequence of the previous one. 
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Lemma 8.4. Let f be a bounded measurable function, with compact support. Sup¬ 
pose that f is positive in an open set U. Then Gm * f is continuous, and super¬ 
harmonic in U. 

Proof. By Lemma |4.8[ and since / is bounded with compact support and the first 
partial derivatives of Gm are locally integrable, it is easy to see that Gm * / is 
continuous. Now, fix a positive function (f G supported in U. Then, by 
Lemma |8.3[ we have 

{Lm{Gm * /), (p) = (-/, (p) < 0 , 

since / and (p are both positive in U. □ 

The next proposition states the convergence of the majorants. 


Proposition 8.5. Let Sn and s the majorants defined on (1.18) and (1.27) respec¬ 
tively. Then 


uniformly on compact sets of . 


By Lemma 8.1 there exists a ball SI C containing supp(m) and supp(M„) for 


all n. Outside SI, we have 


In 


7 = 5 


uniformly on compact sets by Proposition 7.2 We only need to prove convergence 
in S2. 

Let us consider four compact sets Ki, K^, and K^, such that supp(/ 9 ) C Ki, 
SI + Bm +1 C Ki, and 

Kj +Bm+i C Kj+i, for j = 1,2,3. (8.3) 

By our assumption about Ki , Lemma |8.2| and Proposition |5.6| we have that the 
discrete and continuous majorants can be written as 

Sn = inf{/(x); Lm,™/ < 0 on K'j and / > 7 „}, (8.4) 

s = inf{/(x); / is continuous, Lm/ < 0 on Ki and / > 7 }. (8.5) 

Consider the function u}(x) such that 0 < w < 1, uj{x) — 1 for x G K 3 , and 
supp uj C K 3 -\- Bi. 

Define the following function on 7 ^^: 

(pn = -LM,n(s„w). ( 8 . 6 ) 

Lemma 8.6. The sequence of functions {(pn)n is uniformly bounded. 

Proof. Since supp(a;) C + Bi we have, supp((^„) C {K^ + Bm+i) C K 4 . 

li X G K 2 , (pn coincides with —Lm.uSu which is uniformly bounded by 1, as can 
be easily seen as a consequence of the odometer’s equation. 

Consider then x G K 4 \ K 2 ='■ E. In the set E we can replace s„ by 7 „, then in 
(] 8 . 6 |), we have 

(8.7) 

( 8 . 8 ) 

By Lemma 4.3 the last term on the right-hand side of ( | 8 . 8 [ ) converges to Lm( 7 w)(x) 
uniformly on E, once we have that 70 ; is G^. We conclude that LM,n{l<-^){x) is 
uniformly bounded in E. 


We write 


(pn{x) = -LM,ni'JnUj){x) for X G E''. 

LM,n{'lni^){x) = LM,n{{ln “ 7)w)(x) -|- L M ,n{lt^){x) ■ 
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We only need to check that the first term on the right-hand side of (8.8) is also 
uniformly bounded on E. 

Note that 


iM.n((yy -^)w)(a;) = - ■ Pw)(a;) 


l-l 


is uniformly bounded on E, also as a consequence of Lemma |4. 3 
Then, it suffices to prove that 


LM,n{{GM * p — GM,n * Pn)^) 


(8.9) 


is uniformly bounded on E. 
We write 


J^M,n — - n —T-r 

n^-a 


where 


^ f{x + y)-f{x) 




and LM,nf{x) = ^ 5^ 


Fu{y){f{x + y)-f{x)) 


yeB-^UdB-^ 


\y\ 


2+a 


for a function f dehned on -Z^. 

Let us see that LM,n{{GM * p— GM,n * Pn)^) is uniformly bounded on E. Recall 
that for a: G we dehned Pn{x) = p''{x). 

Denote Em '■= E + Bm+i- For x G E'^j we have 


Gm * p{x) - GM,n* p{x) = / {Gm{z)p{x - z) - GM,niy)p{x - y))dz 

1 T2 ■'v° 


= Ii{x)+l 2 {x), 


where 


hix) ■=^ {Gm{z)p{x- z)-GM[y)p(.x-y))dz, 

h{x) ■= \ i‘^M{y)p{x -y)- GM,n{y)Pi^ - y))- 




Since Ki contains the support of the initial distribution p, the sums on the dehnition 
of /i(a;) and hix) are taken over a set which is safely away from the origin. Let us 
say E. 

By symmetry, we write hix) as 



GMiy+z)pix-y-z)+GMiy-z)pix-y+z)-2GMiy)pix-y)dz. 

J-12 


Since the sum above has been taken over a set E that does not intersects a neighbor¬ 
hood of the origin, the function Gm, and consequently the function GMi')pix — •) 
is that means that using second order Taylor’s Theorem we have 
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a^) < ^ 51 sup i:»'^(GM(j/)p(a;-2/) 

^ 1 ^ 1=2 
yeE- 

C 
2^ 


< ^m(£')sup sup \Dl^{GM{y)p{x-y))\ 

yeBl/3|=2 

C 


< 


( 8 . 10 ) 

( 8 . 11 ) 

( 8 . 12 ) 


where p denotes the Lebesgue measure in and C is a positive constant depending 
only on a, which has changed from one line to the other. 

Let us now estimate 12- As noticed before, the sum I 2 (as the sum Ii) is t aken 
over a set E that does not intersects a neighborhood of the origin. By Theorem 1.3 
we have 

\h{x)\ < (8.13) 

n“ 

where G is a positive constant which only depends on a and on the set E. 

We have proved that, for all x G Em 

\Gm * p{x) - GM,n* Pix)\ < (8.14) 


As a consequence of ( 8.14| ), ior x G E we have 
|LM,n((GM * P-Gm.k * p)w)(a;)| < 


G 


^2+a 


E 


\y\ 


2+a 


^ c E 


1 ,, \y\ 


1 

2+a ■ 


y£7P\B. 

which proves that LM,n{{GM * p — GM,n * p)<^){x) is uniformly bounded in E. 

The same argument proves that uniformly bounded on 

E, which concludes the proof. □ 


Lemma 8.7. The function (fn defined in (8.6) satisfies 

|s° -G* ((/>„)°| 

uniformly on K 2 ■ 


0 


Proof. In the beginning of the proof of Lemma 8.6 we observed that supp((/)„) C K^. 
Note that 

Gm * ((/'n)°(x) = 5^ (fniy) [ Gm{x- Z)dz, 


veK, 


and 


s„(a;) = Gm.u * 4>n{x) for x G Kif. 

Hence for x G K 2 

S°{x) - Gm * {(t)n)°{x) = GM,n* (t)n{x-'-) - Gm * {(t>n)°{.x) 

= E ^^^y)^n(.y), 


y&K- 


for a.f^{y) defined in (ini)- 
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By Lemmas | 8 . 6 | and [71^ there exists a positive constant C independent of n and 
X, and a positive function / in such that 

C 


|s„(a;) - Gm * ((/>n) (a;)I < 


,1+a 






< 


C 


— nO!-l 


IK 4 


f{y-x')dy 


< 


C 


nOL — 1 


’{K2+Ka) 


f{y)dy 


< 


c 

tQ :—1 ' 


(8.15) 

(8.16) 
(8.17) 


where the constant C has changed from one line to the other. In (8.151 we used the 
approximation of the integral by Riemann’s sums, in (8.16) we used that x € K 2 
and y € and in (8.17) we used that / is This finishes the proof. □ 

Proof of Proposition \8.5\ As we pointed out before, we only need to prove that s„ 
converges to s uniformly in a ball as the one considered in Proposition 8.1 Recall 
from (8.4) and (|8.5[ ) the expression we are using for s„ and s. 

Recall from (8.3), our assumptions about the sets Kj, for j = 1,2, 3,4. We will 
prove that s„ converges to s uniformly on Ki. 

We write 


s{x) = f s{y)X \(^^^-^)dy, 

JR2 A 

where t] is the standard smooth mollifier 


(8.18) 


il{x) = 



if |a;| < 1 
if |a;| > 1 


normalized so that J ^2 ydx = 1 (see [3] ). Then s is smooth and superharmonic. 

By Proposition |5.6[ s is continuous, and by compactness, s is uniformly contin¬ 


uous on K 2 , so taking A sufficiently small in (8.18), we have |s — s| < e in K 2 . 
By Proposition |4.3[ there exists a positive constant Gi such that the function 


dn{x)=S 


7.2 


we have 7 ^ —>■ 7 uniformly in 7^2 ■ 


satisfies LM,n<ln < 0 in Kf. By Proposition 
Taking n large enough to ensure that < e in K 2 and | 7 „ — 7 ''| < e in Kf, 

we obtain 


qn > s " — e > s " — 2e > 7 " — 2e > 7 „ — 3e 


in Kf. 

Now, the function /„ = max(g„ + 3 e, 7 „) coincides with + 3e in Kf, then /„ 
satisfies LM,nfn < 0 in Kf. It follows that fnP Sn, hence, in Kf, we have 

Sn < + 3e < s'' + 3e < s''' + 4e. 

By the uniform continuity of s on Ki, taking n larger if necessary we have |s—< 
e in Ki , and hence < s + 5e in ATi. 

For the reverse inequality, recall from ( 8 . 6 ) the definition of the function <()„. By 
Lemma |8.7[ we have 

|s° -Gm* ((('„)°| < e 
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in K 2 , and hence 

Gm * > 7° - e > 7 - 2e 

in K 2 for n sufficiently large. Since (/>„ is nonnegative on K 2 , by Lemma |8.4| the 
function Gm * is superharmonic on K 2 , so the function 

Qn = max(GM * {4>n)° + 2 e, 7 ) 

coincides with GM*{(t>n)^+‘^£ in K 2 , then is superharmonic on Ki for sufficiently 
large n. By the definition of s it follows that Qn > s, hence 

> Gm * (</>n)° - e > s - 3e 

on iLi for sufficiently large n. □ 


Proposition 8.8. Let and u respectively be the discrete and continuous odome¬ 
ter function defined in (1.191 and (1.281. Then 


uniformly as n goes to infinity. 



9. Convergence of the final distribution 


Now we are in position to prove the main theorem of this paper, which tells us 
about the weak-* convergence of the final distributions. 


Proof of Theorem \1.S\ Recall from (1.21) that the final distribution of the trun¬ 
cated a-stable divisible sandpile in starting with initial density of mass = p '' 
satisfies the equation 


Let (j) ^ Gf 


Un — Pn ~\~ LM,n‘^n- 
be a test function. We have 


{Pji ; fi) — GMjU * Pn)) ^ 

~ (^n {Gm^u * Pn) 


(9.1) 


Lemma 

4.3 

and 

7.2 

we 


e the convergence (un — Gm.u * p") — t u — Gm * P uniformly on 

compact sets of K^. Since </> has compact support, there exists a ball which contains 


the support of pM.nfi for ri. From (9.1) we have 

lim (v°, (/)) = {u- Gm * P, Lm^) = (p + Lmu, fi). 

n—>-oo 

Let us denote v = p Lmu, we have just proved that converges towards v in 
the weak-* topology. We must show that the distribution 1 / is actually a function. 

By Lemma |8.1| we can also say that there exists a ball U which contains the 
support of the final distribution for all n. Then we have 

(^n = / ^nix)Hx)dx. 

Ju 

Since 0 < < 1 for all n, by Holder’s inequality we have 

\{v,(l))\= lim \{v°,(l>)\< [ |())(a;)|da; </r(C/nsupp((^)) II (() Iloo (9.2) 

n-i-oo Jjj 
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where /i denotes the Lebesgue measure. This shows that p is a continuous linear 
functional, then by the Riesz-Markov Theorem, there exists a Radon measure A 
such that 


= j (t){x)dX{x). 


Moreover, the measure A is non-negative, since the functional v is positive. 

Let A be a measurable set, consider a family of measurable sets such that 

Af: A when e J, 0. Now consider a family {ge)e of infinitely differentiable functions 
such that 0 < < 1 and supp(ge) C A^. In particular, from (9.2 1 , for all e we have 

A(A) < < ^i{A,) II 5 , Iloo . 

Taking limits as e J, 0 we obtain 

A(A1) < ^,{A). 

Then A is absolutely continuous with respect to the Lebesgue measure, therefore 
the linear functional v is indeed a function. □ 


10. Proofs of analytical estimates 

In this section we will proof the propositions and lemmas of Section We will 
often use in our estimates an universal constant C which may change form line to 
line and depends only on the parameter a, except when we specify another possible 
dependency. 


Proof of Lemma \4A\ For R > 0, consider the function Fji in Z^, defined as 

FR{y) = K[y- hv + ^ Fr). 

For a domain A in we denote A"^’'- for its restriction to -Z^, that is A"^’’ 
A n And denote the discrete boundary as 

_ s ^ iZ^; [x — fF x+ ,Ai 2 (t x + ^ 1 ^ n A 7 ^ 0 }. 

^ n ’ I- 2n ’ 2nJ ^ I- 2n ’ 2nJ ' J 


Notice that Fr is equal to 1 in " \ dB 
Let us denote 




]j’ ', and equal to 0 outside B^^ " U dB^". 

Fiy) 


Consider R' > R. A change of variables allows us to write 


\y\" 


k(R!) - k{R) = 


^2-ol 




1 


A^ a 


\y-^ 


2 

|y|“ 


\dz 


We use second order Taylor’s Theorem and obtain 

C 


\k{R:)-k{R)\ < 


-j 2 +q; 


E 


\y\ 


2+a 


< 


< 


c 

n“ 

C 


/ 


R^\Bj 


\y\ 


2+a 


dy 










30 


SUSANA FROMETA AND MILTON JARA 


We have proved that the sequence {k{R)}ii is Cauchy, so it does converge. We can 
ensure then that there exists S > 0 such that 

|fc(ii)| < S for all R. (10-1) 

Although the sequence k{R) does converge, neither the integral nor the sum on its 
definition converge. Then there exists r > 0 so that 


I] ^> 5 - 

.A- 

y&Br-’ 


( 10 . 2 ) 


Taking r as in (10.2) and M bigger that r, by equation (10.1) we have that ki = 
ki{r,M) as in (1.2) is positive. □ 


Proof of Lemma \4-S\ Denote by k := lim„_,,oo kn- In particular, from the proof of 
Lemma |4.1[ we have 

C 

\k-kn\ < —, 


for a positive constant C which only depends on a. 
On the other hand, we have 

Fr.il) 


E 


1 


E 


< 


|y|2+a 


E 

c 


\y\ 


2+0! 


r}Ot 


also for a constant C depending only on a. This concludes the proof. 

Proof of Lemma \4-3[ Let us write 

4c„fc„A„/(a;) c„ Fniy)ifix + y)-fix)) 


□ 


krM,nf (^) — 


E 




\y\ 


2+a 


where A„/(a;) = \ T,y(^xz^,\y\=xifix + y) - fix))- 

Let us suppose x G K, for a compact set K. Since f G there exists a constant 
C, which depends on the first and second derivatives of / in K, such that 

4c„fc„A„/(x) ^ CAfjx) ^ C 


2 — Qi Ti 2 — ct 


Denote 


Since, by Lemma 
prove that 


9xiy) = 


fix + y) + fix-y)-2fix) 


4.2 


|c„ - c| < 


2|?/|2+« 

+ , for C > 0 depending only on a, it suffices to 


Mdy-^ E 


Fniy)9xiy) 




< 


C 


(10.3) 


for a constant C depending on the first and second derivatives of / in AT + Bm- 
In order to prove the inequality (|10.3[) we write the left-hand side as: 


9xiy)dy+- 


E 


g^^y{z)dz, 
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where 5a;,y(z) = gx{y + z) + gx{y - z) -2g^{y). For 5 small, we have \gx{y)\ < |^, 
for X G with C depending on the second partial derivatives of / in K. Then 

f C 

/ gxiy)dy < - 
iBx 


We need to estimate 


I := 


E 


gx,yiz)dz. 


(10.4) 


yeA-i^M^dA-XM' 


We use second order Taylor’s Theorem in every argument of the sum and write 

gx,v{^) = 

|/3|=2 

where Rp{z) = j^{\ — t)Ddg^ y(tz)dt. Since / is twice differentiable, there exists 

a positive constant C which only depends on the first and second partial derivatives 

2 n ’ 2nJ 


of / in iF + and in the parameter a, such that for all z G holds: 


\Rp[z)\ < <^(1 |2+„ + I |2+a)- 


.|5 + zp+“ \y-z\ 

Notice that \z\ < ^ it z G since y € \ {0}, we have (for a bigger 

constant C) 


\M^)\ < 


c 

i\y\ - 


We conclude 


1^1 ^ S 

y 
C 


_ V2)2+o 

2n > 


< 


(|y| 


< ^ 

— 

This finishes our proof. 

Proof of lemma \4.4\ Note that 

LM{GM*(t>){x)= GM{y)LM4>{x - y)dy. 

By selfadjointness of Lm, we have 

Lm{Gm*4'){x) = / LmGm{x - y)(l){y)dy 

JR 2 


□ 


where 


F (y) = 

2(27r)2 


= / Fx{y)dy, 

r j- ^i{9-x-y+z) _|_ ^i(e-x-y+z) _ 2F(^-x-y) 


dOdz. (10.5) 
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Since the integrand in (10.5) is not absolutely integrable, we are not able to apply 


Fubini’s theorem directly. To overcome this difficulty we make the following trick, 

c</'( 2 /) f 


Fx{y) = lim 


t-i-o 2(27r)2 


r i{e-x-y+z) , i{e-x-y-z) _ n i(e-x-y) 

/ dOdz 

( 10 . 6 ) 


Now we use Fubini and integrate first on z. By the definition of ipM we have 


Fx{y) = lirn 




t^o (27r)2 

Hence, using Fubini again 

/ (j){y)LMGM{x - y)dy = lim 


^i(e-x-y)^-t\e\^^g^ 


= lim 


-1 

-1 


t —^0 

-1 f 


'b2 Jr2 




(27r) 

= -'/'(a;)- 

Since ij) S , in particular (j) G L^. This fact justifies the last two lines above. □ 


Proof of Lemma \4-5\ Before starting the proof, a little note about the constants 
may be useful. 

Remark 10.1. Recall that the constant c is defined by c = limc„. Let us denote 
by Ca and Ca, the constant satisfying 

1 — cos(0 • x) 


Ca\0\°' = C 


|2+a 


•dx 


and 


cos{9 ■ x) 


U |a^l“ 

__ 2 

Standard computations shows that Ca = ■ 

If \9\ < 1 we use Taylor’s Theorem and see that 


dx. 


-ipiie) = c 


1 - cos{e ■ y) 


2+a 


3, 2|l/| 


for a smooth function g m. Bi. 

Note that 

c f {0-yy 
‘^Jb, \y\ 


2+a 


dy = 


c\0? 


al\9\^ 


3i |y|' 


dy+ \9\^gii9), 


-dy 


Then for 9 e Bi 


MO) = ^\O\^ + \O\^9ii0)- 


(10.7) 


Since ifiid) = 0 only if 6* = 0, we easily see that there exists a smooth function 52 
in Hi such that 

1 4 


M0) <xj\9\^ 


+ 92{0)- 


( 10 . 8 ) 
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Now let US suppose |0| > 1. Recalling from 10.1 the definition of the constant Cq,, 
we write 

1 — cos (0 • y) 


tpliO) = - C 


\y\ 


2+Oi 


-dy. 


(10.9) 


The function hi dehned by hi{6) := c ^ j^J^+t^^ dy is clearly bounded. We can 
easily see that for \9\ > 1 there exists a bounded function ^ 2 , such that: 

1 1 h2{9) , „ 

Furthermore the function /12 satisfies for i = 1, 2, and A(^^i^) = 

l^ 2 i +2 for bounded functions h 2 ,i and /13 in The proof of this assertion lies in 
the computation of the first to derivatives of ^ |g|a outside Bi. 


We can write Gi as 


1 




where 


Let us denote 


Gi{x) = 


iBi 


cos{9 ■ x) — \ 

V'i(fi') 


d9- 


cos{9 ■ 


Ibi V'i(^) 


-d9. 


( 10 . 11 ) 

( 10 . 12 ) 


hix) = 

f cos( 6 » • a;) - 1 

A. w 

(10.13) 

f2{x) = 

f cos{9 ■ x) 

A: l»l" 

(10.14) 

fsix) = 

j {cos{9 ■ x) — I)g2{9)d9, 

(10.15) 


Jbi 


U{x) = 

f h2{9) cos{9 ■ x) 

A: l«P“ 

(10.16) 


By (10.8 1 and (10.10) we have 


1 


<^ 1 ( 3 ;) = H- f 2 ix) + fsix) + f4ix). 


(10.17) 


We split the estimation of C?i in four steps, consisting in the estimation for |a:| < 1 
and I a; I > 1 of the functions /i, / 2 , /a and / 4 . 

Estimate of fi: Suppose first |a;| < 1, since in (10.13) we are integrating a 
bounded function, fi is also bounded in the unitary ball. That is, there exists a 
constant C > 0 such that 

|/i(x)| < C, for all |x| < 1 . (10.18) 

Now let us assume |x| > 1, we denote i = |fy and perform a change of variables to 
see that: 

/.(X) = I 




\oy 

1 


d9 




cos{9 ■ x) 


d9 


cos{9 ■ x) — 1 

w 


d9. 
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The first term of the sum above can be explicitly computed and give us a logarithmic 
term: 

=-27rlog|a:|. 

Notice that the last term is a constant, let us call it Ci. Let us investigate the 
second term, we call it h{x), that is 

f cos(0 ■ i) 

/ 

We integrate by parts in the following way 
f A(1 — cos{9 ■ £)) 

/ ” 


h(x) = 


w 


-de. 


(10.19) 


h{x) = 


w 

4(1 — cos(0 • x)) 

w 


■dO 


d0- [ cos(6' • x)dS 


-[ i^-COs{0-x))^r^dS, 

where v represents the inward pointing unit normal along \ Bi) (see [3]). 

The integration by parts procedure allows us to conclude that h{x) is bounded. 

We have proved that there exists a explicitly computable constant Ci and a 
bounded function h{x) given by (10.19), such that 

/i(a:) = —27rlog |x| + h(a:) + Cl, for all |a:| > 1. (10.20) 

Estimate of f^'- Recall from |10.T] the definition of the constant £„. Let us write 

cos(0 • x') 


h{x) = 


W 


-d0 

cos(0 • x) 


IB, 


\0Y 


d9. 


Let us investigate the behavior of the function 

cos(0 • x) 


f2ix) := 


iBi 


\0\^ 


dO. 


We easily see that the function /2 is bounded for |a:| < 1. Then there exists a 
constant C > 0 independent of x such that 


f2{x) - 


< C, for all \x\ < 1. 


( 10 . 21 ) 


Let us study its asymptotic behavior for x large. Suppose \x\ > 1, we write 


Mx) = 


A(1 — cos(0 • a:)) 

' 


d0. 


We perform an integration by parts bellow. We again denote by v the inward 
pointing unit normal vector along the surface dBi. 

^ I — cos{0 -x) 1 


hix) = 


a 


Jb^ 


(^a-cosie.xjjL^ds 
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The first two integrals above decay as |^, and the last one decays as |^. Since 
we are considering |a:| > 1 we obtain that there exists a constant C > 0 which only 
depends on a, such that \f 2 {x)\ < 

Since f 2 {x) = — / 2 (t), and since 2 — a < 1, we obtain that there exists a 

positive constant C which depends on a, such that 

\M^)\ < 1^1 - ( 10 . 22 ) 

Estimate of f^: It is easy to see that fsix) is bounded for |ai| < 1, that is, there 
exists a positive constant C depending on a, such that 

|/ 3 (a;)| < C, for all |a;| < 1. (10.23) 

If |a:| > I, we have 

fsix) = [ g 2 i 0 )cos{e ■ x)d0 + C 2 , 

J Bi 

for a constant C 2 . We perform an integration by parts in the first term on the right 
of the above formula, just like we did in the estimation of /2 and conclude that 
there exists a positive constant C depending only on a, such that 

\h{x) - C 2 I < for all 1x1 > 1. (10.24) 

|x| 

Estimate of f^: Since a > 1 and consequently 2a > 2, we easily see that /4 is 
bounded. Then there exists a positive constant C depending on a, such that 


|/ 4 (x)| < C, for all |x| < 1 . 


(10.25) 


Consider \x\ > 1. Recall from the paragraph bellow equation (10.10) that there 
exists bounded functions /13 and /i 2 ,i for i = 1,2 such that A \g\ 2 l = \g\ 2 a +2 and 


d h2{e) _ h2,i{9) 
dOi |6»|2° — |e|2o + i ■ 

We write 


= ITT? 


h 2 { 0 )A{l — cos{0 ■ x)) 


d0. 


Denoting by v the inward pointing unit normal vector along dBf, we integrate by 
parts and write 


fiix) = 


1 


h 3 { 0 ){l - cos{0 ■ x_))^^_ 1 _ 


\0\2a+2 

h2{0) d 


IdBI 


(l-cos(0-x))|^^d5 


\x? JdBl 


(1 — cos(0 • x))dS. 


The first two integrals above decay as and the last one decays as Then 
there exists a positive constant C independent of x such that 

C 


|/4(x)| < 1 ^, for all |x| > 1. 


(10.26) 

We have proved that the functions fi , f^ and f^ are bounded in the unitary ball 
Bi (see (10.18), (10.23) and (10.25)). In (10.21) we proved that the function /2 
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has a special behavior in Bi. By (10.17), we conclude that there exists a positive 


constant C which only depends on a, such that 

n ! \ ^ 

Gi{x) - - 


\2-a 


< c. 


(10.27) 


See Remark 10.1 to justify the constants in the formula above. 

Now, in (10.22), (10.24) and (10.26) we estimate the asymptotic behavior of / 2 , 
fs and / 4 . Since the lowest decay of the three functions is jypTr, which is the decay 
of / 2 , we can say that there exists a constant C > 0 depending on a, such that 

l/2(ai)|,|/3{a:)|,|/4(ai)| < ^ 


|2-ct ■ 


Finally notice that the function fi has a special behavior in Rf (see (10.20)). Then, 
by (10.17) we conclude that, there exists an explicitly computable constant (7, and 
a positive constant C independent of x, such that 


Gi(x) + ^ log |x|- %h{x) - C 


< 


C 


—, for all |x| > 1, 


(10.28) 


where the function h is given by (10.19). 

By (10.11) and since we fixed xq so that |xo| = 1, we conclude the proof. □ 


Proof of Proposition \4.6[ Recall from ( |4.2[ ) and ( |4.3[ ) the definitions of the functions 
g, h and h^. By Proposition |4.5| we have that, for x G Bm^ holds 


(10.29) 


.X “1 X 

^ (27r)2c 


and for x S holds 


G'i(^) = -^ log |x| + ^ logM + hi(^) + 


'M' 


M |x|' 

kol 


'1 "^1 

Since we fixed Xq, so that |xo| = 1, we will always have < 1, then 

,, 2 ^ 2 -a 


, Xq , O' 
- 


I / ^0 

M'~ (27r)2c 


(10.30) 


(10.31) 


By equation ( |4.1[ ) we have 
a2 


Gnix) = 




(2 


TTj^ClXP “ ( 27 r) 2 c M 2 


\2-a 


for X G R 


M 


and 


2 h( — ) 

Gm{x) = - 2 ^ log |x| +6m+ I + hM{x), for xG Bj^ 

'^o'm 


where Sm = — 


3(3?) 


M = - (iFpG + logM - and Hm is given by 


hAiix) = 


cos{9 ■ oj) 




d0. 


for any unitary vector uj. 

In the equations above we use that this can be checked by a 

simple change of variables. □ 
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Proof of Proposition \4.8[ By equation (4.1 1 , it suffices to prove that the first partial 
derivatives of Gi are in 

We recall from (10.11) that Gi{x) = G(x) — Gi(xo), for Gi defined in 


(10.12). Then we only need to prove that the first partial derivatives of Gi are in 

b — 


loc' 


By equation (10.17), it suffices that the first partial derivatives of the functions 
fh / 2 , /a and f^ defined in (10.13), (10.14), (10.15) and (10.16) are in Stan¬ 
dard arguments show this fact. □ 


Proof of Lemma 4-9 Let us write 

Gm*p{x) = / GM{y)p{x - y)dy 


/R2\Bj 


GM{y)p{x - y)dy. 


= / GM{y)p(.x - y)dy 

J Bm 

We consider x big enough so that Bm H B{x, supp (p)) = 0, then we have 
GM*/o(a;)=/ GM{y)p{x - y)dy. 


Now we use Proposition |4.6| to ensure that for |?/| > M, there exists a constant 5m 

(10.32) 


(explicitly computable), and a bounded function Hm in Bf^ given by 

1 


huiy) = 


nr 


cos{0 -u}) 

-^— do, where oj is any unitary vector. 




and a function S{y) converging to zero as y goes to infinity, such that 
Gm^v) = log \y\ + hM{y) + 5 m + S{y). 

Hence, we can write 

GM*p(a;) = / S{y)p{x-y)dy+ {-:,^^og\x-y\+hM{x-y)+5M)p{y)dy. 

Jsunn(p) “ 


Now we differentiate twice, with respect to the variable x, under the integral 
sign. We note that all the second order derivatives of the logarithm function goes 
to zero at infinity, we have to check that the function Lm has the same property, 
we write Hm in polar coordinates: 


1 r 

hM{x) = „ 2 / 


“ cos(rcos/3) 


drdjd. 


This change of variables facilitates the computation of the second order derivatives 
of /iM, so we can easily check that they converge to zero at infinity. 

We have proved that the second order derivatives of Gm * p{x) converge to zero 
as X goes to infinity. 

Notice that the Hessian matrix of the function —is equal to —l-id for all 
X. Then we can choose a ball H big enough, such that the Hessian matrix of the 
obstacle 7 is negative definite in \ H. As a consequence, 7 is concave outside the 
set n. □ 


















38 


SUSANA FROMETA AND MILTON JARA 


Proof of Proposition \4.1^ We consider the concave envelope of the obstacle func¬ 
tion 7 , that is, the least concave function bigger than or equal to 7 . We call this 
function F. 


By Lemma 4.9 there exists a ball O such that the function F coincides with 7 
outside n. We want to prove that F is continuous and superharmonic. 

The concave envelope is always continuous. Since the obstacle 7 is twice differ¬ 
entiable, then F is twice differentiable almost everywhere. 

Fix cc € such that F is twice differentiable in x. Consider a tangent line ix to 
the graph of the function F at the point x. By concavity, £x lies above the function 
F, then T — £x is nonpositive and attains its maximum at x (where it is equal to 
zero). Flence LuffT — £x)[x) < 0. Since the truncated fractional Laplacian of a 
linear function is equal to zero, we obtain that Lm^(.x) < 0 . 

We repeat the same argument for all x such that F is twice differentiable and 
conclude that Lm^{x) < 0 almost everywhere. This implies that for every smooth 
function (j) with compact support and positive, holds 

(r, Lm(I>) < 0. 


We have proved that F is continuous and superharmonic (recall Definition 1.26). 
Flence we can ensure that F > s, where s is the least superharmonic majorant 
defined in (1.271, and since F coincides with 7 in so it does s. Since u = s — 
we finish our proof. □ 


11. Proof of Propositions of Section [6 ] 

In this section we use an universal positive constant C that only depends on a, 
which may change from line to line. 


Proof of Lemma \6.1\ We write 
1 


-IpMiO) - = / 

c Jb, 


1 — cos (0 • w) , 1 

dy - 7 


|y| 


2+a 


E 




K (y )( 1 - cos( 6 l ■ y)) 

\y\2+a 

( 11 . 1 ) 

Since for 6 G Bj_ and y G Bm we have |(d • y)| < 1, by Taylor’s Theorem there 
exists a smooth function g defined in the unitary ball that satisfies 


1 - cos( 6 » • y) = + ( 6 ) . y)‘^g{e ■ y). 


( 11 . 2 ) 


An easy computation shows that 

1 f {e-yf 


\2+a 


Analogously 


E 

yaB-M^OBM 


Bm Ijfl 

Fn{y){0-yf 




\y\ 


2+a 


\yY 




Fnijj) 


4 ^ lyh 


Then using the definition of kn (recall from (1.12)), in (11.1) we have 


-f’MiO) - 1pM,nid) = 
c n 


kn\d\" 


+ Ai. 


(11.3) 
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where 


:= 


[9 ■ yYg{e ■ y) 




\y\ 


2+a 




Fniy)iO ■ y)^9{0 ■ y) 


\y\ 


2+a 


(11.4) 


yeB'^UdB'^ 

Let us denote fy^z) = ■ We write 

{e-z)^g{e-z) 1 


+ = 


|2+a -d^+2 


E 




J^12 

2n' 2ni 


\v\^ 

fy^e{z)dz. (11-5) 


We use second order Taylor’s Theorem with remainder on( 2 ;), since 2 ; G [“ + ? 
we have 

fvAz) = E 

|/3|=2 

where ^(z) = fy^ff(tz')dt. We easily see that there exists a 

constant C which depends on a such that 

1 1 


\Rly/z)\<C\9\'^{ 


\y + z\°‘ \y-z\ 


:r 


Since s e ^1^ we have |c| < ^. then 

i«j (c)|< 

■ (l!/l - #)■ 

for a constant C which only depends on a. Then 
1 


E 


yGAl’-^USAT" ^ 2„.2„ 


l/.WMv < tW 

i_ j_i2 n 


E 


yeAT'' uaAT'V 2n ) 


< Wy- 

- .^ 4 -a 


- ,M ^ ,M 

1 

_ 

2 cl 


< 




9 ’ 


where the constant C have changed from line to line. 


The first term in the right hand side of (11.5) is easily seen to be < +=W for a 
constant C depending only on a. Since 4 — a > 2 we conclude 

ci^i^ 


l+l< 


n2 • 


We have just proved that there exists a positive constant C which depends on a 
such that 


Em(6') - ipM,n{d) - 


kn\8\^ / c\9\ 


< 


which together with (6.4) concludes the proof of the lemma. 


□ 


Proof of Lemma \67^ We consider the function in -Z^ defined as 
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We split the difference — i^M,n{0) in two parts and estimate each one 

separately. We write 


-iPm{ 9) — i’M,n{d) = Bi+ B2. 


where 

Bi ^ 

and 

B2 ■■= f 

J A 1 
W 


l-cos(6»-y) 1 ^ 

.|2+a J.2 2-^ 


F{y){l - cos(6» • y)) 




\y? 


m m 


\v\ 


2+a 


( 11 . 6 ) 


(11.7) 


1 — cos(0 • y) , 1 

-dy —T 


\y\ 


2+a 


E 


{Fn{y) - F{y)){l - cos(6» • y)) 


yeAY' 

W“ FT' 


\y\ 


2+a 


recalling that A i = Bm \ B i . 

|e |' |e| 

Estimate of Bi : For y £ B i , we can write 
"ft 

1 - cos(6» • y) = + {9 ■ y)^g{e ■ y) 

for a smooth function g defined in the unitary ball. Hence, we write 

Bi = Bi^i + Bi^2 

where 

\o?(f %)' 


( 11 . 8 ) 


^1,1 = 


1 \y\ 


E 




—' Ivb 

TF T«T 


(11.9) 

( 11 . 10 ) 


and 


7?1,2 = 


{e-yYg{e-y) 


TF 


\y\ 


2+q 


E 


F{y){e-yfg{e-y) 


TF TF 


\y\ 


2+q: 


( 11 . 11 ) 


We estimate first Hip. Note from the definition of /c„ (see (1.12)) that 

1 

|y| 


fc„|0p |0|+/• 1 1 

-SiA =-F-r --r / TTydy--y 


E 


Fn{y) - F{y) 


W'" yeeAp 

TeT’" T«T' 

Let us denote fy{z) := + |^+p - |^)- Then 


\yY 


( 11 . 12 ) 


„ knW W ^ 

yGA":'- uOA"’" 
T«T’" TeT’ 


J_ J_12 

27x ’ 2n J 


fy{z)dz. 


We use Taylor’s Theorem and bound fy{z) by its second order derivatives just like 
we did in the proof of Lemma |6.1[ We obtain that there exists a constant C which 
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only depends on a such that 


/ ify(^)idz < Cjzj^ sup ( i ,^, 2 +a + I- 


< 






where the constant C have changed from one line to the other. 


\el 

4 


E 


ye A":'-' USA"-;-' 


2n ’ 2 ti J 


[- 2 k’2^]^ we 

have 

\fy{z)\dz < 

c\e\ 

< 

cm 


•n? 

< 

cm 


2n ■ 


E 


1 


r, r, ■■ (\V\ - ^12 + 0 

VeA^:- USA’',’- 2nl 

TST'“ FT'" 


1 


\y\ 


2+a 


dy 


Ter 


2 + Q: 




We have proved that there exists a constant C which only depends on a, such that 

fc„|d|2 C|0|2+“ 


Bia- 


',2 — 0; 


< 


(11.13) 


Now let us estimate Bi, 2 - We proceed in the same way as in the estimation of 
Bia. Using second order Taylor’s Theorem we obtain that there exists a constant 
C depending only on a, such that 


\BiM < 

< 


C|6»|^ 


c\e\^ 


E 


yeB"'{''UdB"'{'‘ 

T«T Wi 

TF 


1 


(M- 


V^ \OL 
2n X 


where the constant C have changed form one line to the other. Therefore 


1 ^ 1 , 2 ! < 


c\e\ 


2+a 


Putting (11.13) and (11.14) in (11.9), we obtain 

C|6i|2+“ 


\B. - 


kn\OY I ^ 

n2-a ' — 


(11.14) 


(11.15) 


Estimate of B 2 : We repeat the same argument used in the estimate of Bi. That 
is, first we write the difference between the integral and the sum, as a sum of 
integrals over small squares. Then we symmetrize the arguments in the integrals 
and use second order Taylor’s Theorem with remainder. We conclude that there 
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exists a constant C which only depends on a, such that 

C\e\^ 1 


\B, 


2 | < 


< 


E 


^ -L- A/f -L- 


_ V2\2+n 

2n > 


c\e\^ 


\y\ 


2+a 


lei 

dy, 


m 


where the constant C have changed from one line to the other. We conclude 

C'|6l|2+“ 


\B2\ < 


The estimations of Bi and B 2 together give us 

= 

c “ 


2+a 


(11.16) 


(11.17) 


where r„ is an uniformly bounded sequence of functions defined in [—nir, rnr] ^ \ . 

Finally, by (6.4) we conclude the proof of the first part of Lemma 6.2 The proof 
of (i) and (ii) follows that same idea used in the first part of the proof. □ 

Proof of Lemma \6.3\ Let us introduce the notation 

f cos{8 ■ x) — 1 


hix) 

/i,«(a;) 

/2(+ 

/2,n(a;) 

/3(+ 


cos{9 • a;) — 1 


-de 


d0 


i’M,nid) 

cos{6 ■ x) 
1pM(d) 
cos{9 ■ x) 


[—nTT,n7r]^\B i 

'M 


[—n-K ,n7T]^\B 1 i’M,n{6) 

AT 

j' cos{9 ■ x) 

'ipM{9) 


d9 

d9 


d9. 


— nTT.nTT ' 


Then 


— GM{x)-GM,n{^) = { — fl{x)-fl,n{^)) + { — f2{x)-f2,n{x))+ — Mx)- (11.18) 


We split the estimation of fhCivi — GM,n in three different steps, consisting in the 
estimation of +/i(a:) - fiAx), +/ 2 (a;) - / 2 ,n(a;) and ^faix) separately. 

Estimate of +/i(a:)— /i,n(a;): ByLemmaf 
G which only depends on a, such that 


6.1 


we have that there exists a constant 


— fi{x) - fi.Ux) 


^ ^ [ (l-cos(6>-a;))|6>|^ 


d9. 


Since ipM,n and ipM are both of the order of jdp in i?j_, we conclude 


— fi{x) - fi,n{x) 
Cn 


c 

< 


(11.19) 






















SCALING LIMIT FOR A LONG-RANGE DIVISIBLE SANDPILE 


43 


Estimate of ■^f 2 {x) — f 2 ,n{x)'. In the region [—n7r,n7r]^ \ the functions 
and tpM are both of the order \9\°‘. 

Denote 0„(0) := ^ ~ dM^~T6)' Lemma 


6.2 


we have 


(l>n[d) = 


^mo\ 


2+a 


■nfCn 1pM{ff)'tpM,n{d)' 


for a uniformly bounded sequence of function r„ defined in [—n7r,n7r]^ \ i?j_. We 
have 


f2{x) - f2,n{^) = J 


[—n7r,n7r]2\B 

M 

-1 f 


1^1 d[—n7r,n7r]^\R i 


COs(d • x)4>n{d)d9 

A(cos(0 • x))(fn{Q)d9. 


Integrating by parts, we obtain 


— f 2 {x) - f 2 ,n{x) = [ C0S{9 ■ x)A(j)ni9)d9 

1^1 J [—n77’,n7r]‘^\B i 
W 


]—fT / cos{9 ■ x) — (j)n{d)dS 


roT / (l)n{9)— cos{9 ■ x)dS, 

1^1 J d{[—n7T^n7T]‘^\B j^) 

M 


where u represents the inward pointing unit norma l vec tor along the surface 
9([—RTF, n7r]^\i?^). By items (i) and {ii) of Lemma 6.2 we have |A0„(0)| < 
then 




I 


' [—n7r,n7r]2\S i 


cos{9 ■ x)A(j)n{9)d9 


< 


< 


< 


C 


1 


|2 I, „ rB\9\°‘ 


d9 


1^1 J \—mT,n'K\2\B 1 


c 


— nTT.nTT • 


I 

w 


-d9 


C 


( 11 . 20 ) 


Now let us estimate the boundary terms. Since 'ipM,n is periodic with period 
[—n7r,n7r]^ we have 


[ cos{9-x)-^ ^ dS = [ ^ IW 

/a[-n7r,n7r]2 012 1pM,n{d) Jd[-n-n,n-^Y VM,n[d) OV 


On the other hand it is easy to see that 


1 


/ cos{9 ■ x)—(t)n{9)dS 

IdB 1 012 


< 


C 


cos{9 ■ x)dS = 0. 

( 11 . 21 ) 

( 11 . 22 ) 


/ 4>ni9)-^ cos{9 ■ x)dS 

IdB 1 012 

-M 


< 


c 


(11.23) 


and 
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Define 

Tn{x) := 


cos{9-x) 


d 


-dS- 


d 


Cn|^| J d[—n-K dvipMiO) 'tpM{d)du 

Then, from (11.20), (11.21), (11.22) and ( 11.23| ), we conclude 

—/ 2 (a;) - f 2 ,n{x) - Tn{x) 

Cn 

Estimate of ^faix): Note that 


c 

' 1 1 \ 

< — 


“ n“ 

a; p T ) 


cos{9-x)dS. (11.24) 


(11.25) 


f/3(^) = f 

Cn Cn . 


cos(0 • x) 


d9 


'R2\[-r 


f’Mid) 

A(cos(0 • x)) 


/B2\[-r 


IpAlid) 


d9 


Performing again an integration by parts, we obtain 


Mx) = 


CnX 


f cos{9-x)A( ^ )d9-Tn{x), 

7R2\[-n7r,n7r]2 


where T„ was defined in (11.24). 
vVdlMP - W 

1 


have that 


on \ i3j_ an in particular on \ [—n7r,n7r]^, we 


1 


n7r,n7r]^ 




Then 


— h{x) + Tn{x) 

CtJ. 


c 

< 


< 


c 


(11.26) 


By (11.18), (11.19), (11.25) and (11.26) we conclude the proof of Lemma |6.3| □ 

Acknowledgments. The authors would like to thank to Charles Smart, Luis 
Caffarelli and Luis Silvestre for fruitful discussions. 
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